Dynamics of a two-species Bose-Einstein condensate in a double well 
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We study the dynamics of a two-species Bose-Einstein condensate in a double well. Such a system 
is characterized by the intra-species and inter-species s-wave scattering as well as the Josephson 
tunnelling between the two wells and the population transfer between the two species. We investigate 
the dynamics for some interesting regimes and present numerical results to support our conclusions. 
In the case of vanishing intra-species scattering lengths and a weak inter-species scattering length, 
we find collapses and revivals in the population dynamics. A possible experimental implementation 
of our proposal is briefly discussed. 



I. INTRODUCTION 

Soon after the realization of scalar Bose-Einstein con- 
densates (BECs) in laboratories, there has been great 
interest in exploring two species BECs both experimen- 
[ la B, M, M , Hol. j n i Il2| and theoretically 
161 . 117L [H, [ljj, [2O, UM ■ Such a two species 
BEC provides an ideal platform for the study of more 
intriguing phenomena. It can help us understand some 
basic problems in ultralow temperature physics like in- 
terpenetrating 3 He- 4 He mixtures. The two species can 
be two alkali like 23 Na- 87 Rb or two isotopes like 85 Rb- 
87 Rb. It can also be two hyperfine states of the same 
alkali. Different from the case of a scalar condensate, 
we need three scattering lengths to characterize the s- 
wave scattering between alike and unlike atoms, which 
is referred to as intra-species and inter-species scatter- 
ing, respectively. The interplay between the inter-species 
and intra-species scattering has a direct consequence on 
the pro per ties of the condensates, e.g. the density pro- 
file [I3LTT5I] and collective excitations [l6[. For example, 
as shown in a pioneering work by Ho and Shenoy [13j |. 
the phase of trapped condensates can smoothly evolve 
from interpenetrating to separate by changing the atom 
numbers and/or the scattering lengths. 

In addition to studying the stationary properties of 
BECs, the coherent dynamics of a two species BEC can 
be investigated by propagating coupled Gross-Pitaevskii 
equations. Previous work has shown various interesting 
dynamics for two-species BECs [H [H, US HI ■ For this 
paper, we focus on the case of the two species being two 
hyperfine states of the same alkali. One can manipulate 
such a two species BEC either by turning on the pop- 
ulation transfer between the two internal states using a 
coupling microwave field [ISIIEIIzljlHilIflllO] or by plac- 
ing it in a double well to allow for Josephson tunnelling 
[I]], H HI [H, HI- For the latter case, the double well 
is fundamental to the study of the Josephson effect in 
BECs. Smerzi and co-workers have shown that a scalar 
BEC trapped in a double well exhibits rich dynamics 
such as Rabi oscillations and macroscopic quantum self- 
trapping [36l . l37l . l38j . It has also been shown that there 
exists non-Abelian Josephson effect between two F = 2 
spinor BECs in double optical traps (39|. Therefore it is 
of great interest to study the behavior of a two species 



BEC in a double well. Previous work on a two species 
BEC has focused either on the population transfer or on 
Josephson tunnelling. To the best of our knowledge, the 
combination of both of them has not been addressed in 
the literature. It is the purpose of this paper to discuss 
some interesting phenomena involving both population 
transfer and Josephson tunnelling. 

This paper is organized as follows. In Sec. II, we first 
introduce our model Hamiltonian and various parame- 
ters. In Sec. Ill, we then discuss some interesting regimes 
and give our numerical results. We summarize our find- 
ings in Sec. IV. Finally, an appendix is given in Sec. 
V. 



II. MODEL 

In this paper, we consider a two species Bose-Einstein 
condensate in a double well. We focus on the case where 
the two species are for two different hyperfine spin states, 
e.g. 87 Rb condensates in |1) = \F = l,mp = —1) and 
1 2) = \F = 2, nip = I), and they are subjected to the 
same trapping potential. We adopt the well-known two- 



mode approximation, i.e. ^ a j(r) 



,^a(r). The 



Greek letter a = L(R) denotes the left (right) well and 
the Roman letter j = 1(2) labels the two species. V'a(r) is 
the ground state solution of the a well which is indepen- 
dent of the species. a a j is the annihilation operator for 
species i in the a well, satisfying [a Q j,ot t ] = 5 a ,p5j,k- 
Here 8 p ^ q is the Kronecker delta. In the second quantiza- 
tion formalism, the Hamiltonian is (h = 1) 

H = ^ '■>"n.j"n.j I "].. ," K.J I >• I. ..>•)..,) 

a=L,R;j = l,2 j 

" («) „(«) 

ES11 t t "22 t t 

a L 

+ 9l2 )a l, l a i,2 a a,iaa., 2 + Sl( a l,l a a,2 + aa,l a a,a) 

e a is the single particle energy in the a well. 2J = 
fdrf L (r)(T+V(r))^ R (r)+fdv^ L (r)(T+V(r))r R (r) is 
the intra-species Josephson tunnelling rate between the 
left and right well which is referred to as tunnelling in 
the following. Here T and V are the kinetic and poten- 
tial (double well) operator of a single particle, respec- 
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tively. J can be tuned by changing the shape of the 
double well, g^j = (Airh 2 'aij j 'to) Jdr|^ Q (r)| 4 is the in- 
teraction strength between species i and j in the a well. 

is the s-wave scattering length between species i and 
j. In this study, we only consider the case of repulsive 
interactions, i.e. > 0. m is the atom mass. 2f2 is the 
Rabi frequency for the inter-species population transfer 
which is referred to as coupling in the following. It is pro- 
portional to the amplitude of the microwave field which 
resonantly couples the two internal states. For simplicity, 
we assume J > and fi > 0. In addition, we consider a 
symmetric double well in this paper, thus cl — e R and 
V>l(i") = V'fl( r )- Correspondingly, g£' = g\^' and the su- 
perscript is then dropped afterwards. A schematic view 
of our model system is shown in Fig. [TJ We note that 
the Hamiltonian of Eq. (TTJ) resembles the cross phase 
modulation between two modes of light in a nonlinear 
Kerr medium. The terms involving Q and J are analo- 
gous to the Hamiltonian of a beam splitter. In the point 
view of quantum optics, our Hamiltonian describes the 
combination of nonlinear beam splitters. 
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FIG. 1: (Color online) A schematic view of our model system 
where the solid lines represent the four modes and the dashed 
lines represents the coupling (red) and tunnelling (blue). The 
s-wave scattering terms are not shown in this figure. 

We assume that initially all N atoms are of species 
1 and localized in the left well, i.e. \ip(t — 0)) oc 
(at 1 ) JV |vac). The coupling f2 and J are turned on at 
t = and we investigate the subsequent dynamics. 



III. VARIOUS REGIMES 

We categorize our model into three regimes: the strong 
coupling/tunnelling regime (fi, J) » 9ij(N — 1), the 
strong interaction regime g^ \N— 1) 3> J), and the in- 
termediate regime where the terms of gij(N — 1), O, and 
J have equally important contributions to the dynam- 



ics. In the following, we will discuss these three regimes 
separately. 

We first consider the strong coupling/ tunnelling regime 
so the collisional interaction terms are neglected. The 
Hamiltonian is 



= ^((2^ -^£,2 + CLr i a R,2 + a L, i<4, 2 + a R,i a R 2) 



+ J( a L.l a R,l + a L.2 a R? + a L,l a R,l + a L,2a Rt2 )- 

(2) 

It is easy to verify that the terms in the bracket involving 
f2 commute with those involving J . Hq can be writ- 
ten in a more compact form = • M • a where 
a = (<Xl,i, o-l,2, aR.ii Ofl,2) T - The matrix M depends only 
on f2 and J. This Hamiltonian is quadratic, so it can be 
simplified in a standard manner by diagonalizing M. Ex- 
plicitly, by the following transformation, 

l r 

;(<i£,i + «l : 2 + a R i 



h 
h 
91 

92 



(ffli,l - «L,2 - OR,l 

{a L ,i + a Li2 - Or,i 
(ffli.i - «i,2 + a R i 



a R,2), 
OR,2), 
OR,2), 



(3) 



[fj, fl] = Sjk and [fflj,o|] = Sjk with all other commuta- 
tors being zero. The conservation of atom numbers leads 
to fth + Ah + 9i9i + 92.92 = N . The Hamiltonian can 
be rewritten in the form of new operators 

H c = (n + J){.flh - flf 2 ) + (fi - J)(g\ gi - g\g 2 ), (4) 

where constant terms have been dropped out. We note 
that the net effect of coupling and tunnelling can be ei- 
ther constructive (fi + J) where the coupling and tun- 
nelling are in phase or destructive (fi — J) where they 
are out of phase, analogous to the center of mass and the 
relative motion of two particles, respectively. The corre- 



0)> «[/} + /. 



-9l 



sponding initial condition is \ip(t 
I vac). 

The wave function at any time t can be obtained with 
the help of the following identity 



-iH c tl A 



ft + 91 + 91)* 
_ e -i(n+j)t yrt _|_ e *(fi+j)t jt + e -i(n-j)t g \ + e »(n- 

Therefore, 

1 



iH c t 



J)t 9l 



\m) 



-i(Q- 



2 JV 
-J)t 



-i(Q+J)t 



fl 



= *(0+ J)t j\ 



9i 



N 



vac) 



Substituting a' a ■ back into the above expression, we ob- 



tain the normalized wave function 



m)) 



=i,H;j=l,2 



N 



3 



S L ,2(t) 
SR,l(t) 

Various quantities can be calculated directly from the 
wave function. For instance, the first order coherence 
between modes a a ,j and ap t k is found to be 

= S* a Jt)S Mt)- (6) 

Especially, the fraction of species j in the a well is given 
by 

n a At)=& a<j (t) = \S a , j (t)\ 2 . (7) 

Higher order coherence terms can be derived straightfor- 
wardly with similar patterns and will not be presented 
here. We note that at any time t > 0, there is a balanc- 
ing condition riLi/n^ — nm/rifa for the given initial 
state. 

Our conclusion also applies to the case with temporal 
modulations of 12 and J if we make the replacement 

e -iHot e -i J* H c {t')dt> ^ ^ 

The time ordered exponential is not needed here since 
the hamiltonian commutes at any time in this represen- 
tation. Therefore, the effect of time dependence is de- 
termined from two pulse areas, i.e. J Q (Cl(t') + J(t'))dt' 
and (fl(t') — J{t'))dt' . If we assume a form of il(t') = 
12oe~* l T and J = Joe - * / r , then the two pulse areas are 
(12 + J )t(1 - e"*/ r ) and (fi - Jo)t(1 - e"'/ 7 ") for the 
case of in phase and out of phase, respectively. For a long 
enough time t> t, the exponential terms in the above 
pulse areas can be neglected and the tunnelling dynamics 
are suppressed. The wave function takes the same form 
as Eq. [5] except for the replacements 12 — > 12o, J — > Jo, 
and t — > t. A numerical example is shown in Fig. [5] Un- 
less otherwise specified, the energy and time scale are ojq 
and l/u>o, respectively, where ojq is the angular frequency 
obtained by approximating one well as harmonic. 

In the opposite regime where the nonlinear interac- 
tion dominates, the coupling between the two species is 
still on resonance. However, the tunnelling between the 
left and right well is frozen due to the large detuning. 
For instance, the energy difference between the two Fock 
states \N) Ll \Q) L2 \0) R1 \Q) R2 and |iV-l)Li|0) L2 |l) m |0) fl2 
is cx gn(N — 1) ^> J for large N. So the dynamics are 
restricted to the left well which is essentially a two-level 
system. A numerical example is shown in Fig. [3] where 
we can clearly see that the dynamics are restricted to the 





— ■ L1 


1 


L2 


- i 


R-| 


t 
i 


R2 


'. A 
















$ \ 
















i >.■ / ^ " 













1 2 3 4 5 

Time 



FIG. 2: (Color online) Fraction of each component as a func- 
tion of time. 12 = 3.2, J = 1.3, r = 7r/4, gn = 322 = 312 = 0. 
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FIG. 3: (Color online) Fraction of each component as a func- 
tion of time. 12 = 3.2, J = 1.3, r = 7r/4, gn = 322 = 312 = 35. 

two species in the left well and the tunnelling to the right 
well is suppressed. 

In the intermediate regime, we cannot obtain solutions 
in closed form in general. In the following, we will con- 
sider a specific case where we can obtain some analytical 
results. We assume gn = #22 = and the term involv- 
ing gi2 is a perturbation to the coupling/tunnelling. We 
denote 312 = g in the following. Motivated by the strong 
coupling/tunnelling case, we adopt the following trans- 
formation 



Oi,l = 


\ 


e iw +*/i + e 


-iuj+t £ 
J2 


+ e ZUJ - 


% 


+ e~ lu} - 




«L,2 = 


\ 


e iw+t /i - e 


-ILJ + t £ 

J2 


+ e ZUJ - 


% 


- e~ lul - 




OR,l = 


\ 


e iw+t /i - e 


-ILJ + t £ 

J2 


- e iul - 


l 9i 


+ e~ lu} - 




OR,2 = 


1 

2 


[e lu,+t fi + e 


-ILU-l-t £ 

J 2 


- e iuJ 


-'91 


- e~ luJ 


-'92], (9) 



where u± = 2(12 ± J), fj and g^ obey the same com- 
mutation conditions as in the strong coupling/tunnelling 



cos (12 + J)t + cos(12 - J)t 



cos(12 + J)t — cos(12 — J)t 



0.8 



. sin(12 + J)t + sin(12 - J)t 



0.6 



.sin(12 + J)t - sin(12 - J)t 



(5) 
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case. Substituting a ai j into the Hamiltonian ([T]), we then 
invoke the rotating wave approximation (RWA) to elim- 
inate the fast oscillating terms and only keep resonant 
terms. Explicitly, RWA requires g < (SI, J, |Q- J|, |2fi- 
J|, |S1 — 2J|). The new Hamiltonian becomes 

Hr = § (flflhh + Mhf2 + g\g\gi9i 



(10) 



+92929292 + 4/ 1 t <7l/i.gi + 4/ 2 f 5 |/ 2fl r 2 
+ (Jl + J)(/ 1 t / 1 -/t/ 2 ) 
+{n-J){g{ gi -glg2). 



This Hamiltonian is separable in two groups of modes 
(fi,g%) and {f-2,92)- Each group of modes, fj and gj, are 
coupled by the mean field interaction and the detuning 
between them is 2J. The above Hamiltonian is already 
diagonal in the basis of fj and gj, so we label the Fock 
eigenstate as \mnpq) in the basis order of /i, / 2 , gi, gi- 
The corresponding energy is E mnpq = E m , v + E n>q with 



E,, 



E n 



(m 2 — m + p 2 — p + Amp) 



+(Q + J)m + (SI - J)p, 

0/_2 _ , 2 



n I 2 

8 (n 



//, + g~ — g + 4ng) 
-(0+ J)n- (SI- J)g. 



Note that as |vac) = |vac), the wave function can be 
computed in the new basis as 



\m) 

„-iH R t 



2 N VW. 

lHRt w E 6 



Ul + fl+g\ + gl) N \^c) 



2^V / j "N,m+n+p+q 
mnpq 



yjm\n\p\q\ 



\mnpq) 



~ 2N E 



n+p+<J 



mnpq 



y/m\n\p\q\ 



" lpqt \mnpq). 

(11) 



The fraction of each species n a j can be computed with 
the help of the following results where the average is 
taken with respect to the wave function \ip(t)}: 



(flh) 

(flh) 

(flgi) 
(AW 

(fW) 



{flh) = (g\gi) = (9I92) = j, 

2i(n + j) t N ( cos(tg/4) + cos(tg/2) 
4 ^ 2 

iV-l 



iV-l 



„2iJt - 



N f l + cos(tg/A) 



nt N f cos(tg/4) + cos{tg/2) 



-wnt-^ /cos(tff/4) + cos(*5/2) 



AT-l 



JV-l 



(/ 2 f 92) 
(Sl02> 



-2iJ 



Jt iV / 1 + cos(i 5 /4) 



JV-l 



4 V 2 
.2i(n-j)t^ ^ cos(t. 9 /4)+cosftg/2) \ 



iV-l 



fl2) 



In the Appendix, we give a derivation of (flh)- Other 
average values can be derived along similar lines. 
After some algebra, we obtain 



1 



n L ,i 



1 + cos(2 Jt) 



1 + cos 



(*g/4) \ 



JV-l 



+ cos(2Sli)(l + cos(2Ji)) 



cos(tg/4) + cos(tg/2) 



jv-r 



(13) 



n L .2 



1 + cos(2 Jt) 



1 + cos 



(*g/4) \ 



JV-l 



cos(2Sli)(l + cos(2Jt)) 



cos(tg/4) + cos(tg/2)\ 
2 J 



JV-l" 



(14) 



1 - cos(2 Jt) 



1 + cos(<5/4) 



JV-l 



-cos(2fii)(l - cos(2 Jt)) 



cos(tg/4) + cos(tg/2) 



JV-l" 



(15) 



n R ,2 = - 



1 - cos(2 Jt) 



l+cos((t g /4)) \ 



JV-l 



9 (20t)(l - cos(2Jt)) ^os(t. 9 /4) + cos(t g /2)^ 



JV-l" 



(16) 



Therefore, in the case of the vanishing intra-species scat- 
tering, n a j will be subjected to collapses and revivals 
(CR). CR is a quantum mechanical effect which is well 
known in quantum optics. It is also found for a scalar 
BEC [40|as well as for a two-species condensate in a sin- 
gle well 27| . The nonlinearity g determines the envelope 
of the revivals as well as the time separation between the 
adjacent collapse and revival. It is easy to show that the 
width of the revival in the time domain is cx l/(gyN — 1) 
for N > 1 and the separation between a neighboring CR 
is 8Tr/g. The coupling SI and tunnelling J determine the 
detailed structures of the oscillation inside the revival 
envelopes. Here we want to point out a major difference 
between our study and that of a two species BEC in a 
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single well. In our case, CR is not observed for general 
scattering lengths a^-. When the intra-species scatter- 
ing lengths are not negligible, the population dynamics 
usually display quite complicated temporal patterns. We 
note that, for large N, since n a j are all small except for 
gt/4= 2mt, we have 



1 + cos(i£/-/4) 



,-sin 2 (t ff /8) 



cos(tg/4) + cos(tg/2) 



-5sin 2 (tg/8) 



(17) 



So up to some factors, our expressions take the same 
exponential form as those in Ref. [27|, e.g. Eq. (31). 
This is expected because the relative difference between 
a coherent state and a Fock state is small for large (aver- 
age) atom numbers. We also note that, for N — 1, n a j 
reduces to the results of the strong coupling/tunnelling 
regime. 

The population difference 
wells/species takes a relatively 
is shown below 



between the two 
simpler form which 



D 



LR 



(n L ,i + n L ^) - (n Ri i + n R ,i) 

N-l 



cos(2Ji) 



1 



3s(W4) V 



Da 



cos(2n*) 



cos(tg/4) + cos(tg/2)\ 
2 J 



N-l 



.(18) 



Again, both of them show CR which are modulated by 
a sinusoidal oscillation with individual frequency 2J or 
20, i.e. J (f2) has no effect on £> 12 {Dl R ). This is not 
a surprising result even in the presence of nonlinearity 
since we are dealing with a symmetric double well. 

A numerical example of ri^i as function of time is 
shown in Fig. |4] where the (red) solid curve is for re- 
sult of Eq. (TTB"|) and the (blue) dotted curve is for exact 
numerical simulations with the full Hamiltonian of Eq. 
(HJ). We can see that RWA does capture the essential 
features of the dynamics. It reproduces the results very 
accurately for the initial time. However, it is only quali- 
tatively accurate for long times, e.g. RWA overestimates 
the peaks of population oscillation for the first revival. 
Both RWA and exact numerical simulations predict the 
CR phenomenon. Therefore, CR for a two species BEC 
is not an artifact of RWA, but rather a consequence of 
phase coherence under nonlinear interaction. 

The general case can be attacked numerically. In Fig. 
[51 we show our numerical results for til,i as function of 
time. The parameters are chosen so as to fall into the 
intermediate regime. For comparison, we also show the 
results of the strong coupling-tunneling case. We can 
see that, with appropriate nonlinearity, the dynamics of 
n^i exhibit completely different and more complicated 
patterns than the strong coupling-tunneling case. The 
nonlinearity not only reduces the oscillation amplitude 
but also destroys the periodicity (4l| . This reduction in 
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FIG. 4: (Color online) riL,i as a function of time. The red 
solid curve is for result Eq. 1)130 , The blue dotted curve is 
for exact numerical simulations with the full Hamiltonian of 
Eq. (HJ. Q = 2.4, J = 0.6, g n = g 22 = 0,512 = g = 0.2, and 
N — 6. Time is in units of 11 q. 



the amplitude is due to the effective dampingwhich is 
common feature in interacting condensates [42j, |43j . 




FIG. 5: ul,i as a function of time using exact numerical sim- 
ulations with the full Hamiltonian of Eq. (JTJ. Q = 1, J = 1.7, 
N = 6. The red solid curve is for gll — g22 — 0.45 and 
gl2 = 0.3. The blue dotted curve is for the noninteracting 
case gll = g22 = gl2 = 0. 

Now we briefly discuss a possible experimental imple- 
mentation of our proposal. Our model can be conve- 
niently realized in an optical super-lattice by superim- 
posing two-color lasers EH • Each lattice site consists 
of a double well potential. J is determined from the shape 
of the double well which can be tuned by engineering the 
ratio of the two laser intensities. The s-wave scattering 
gij can be tuned by using the well established Fcshbach 
resonance. f2 can be tuned by changing the intensities of 
two microwave/RF fields both of which are coupled to the 
third intermediate level with a large detuning, which can 
be adiabatically eliminated. The relative strength in our 
model Hamiltonian can thus be tuned over a wide range. 
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For the CR as shown in Fig. [4] we give an order of magni- 
tude estimate of typical experimental parameters here. If 
we approximate each lattice site as a harmonic trap with 
trapping frequency lu ~ 10kHz. The intra-species inter- 
action strength g can be tuned to be ~ 100Hz with the 
corresponding a±2 — 25as (a-B is the Bohr radius). The 
Josephson tunnelling rate J and the coupling strength f2 
is 300Hz and 1.2kHz, respectively, both of which can be 
implemented using the state-of-the-art technology. 



IV. CONCLUSION 

To conclude, we have studied the dynamics of a two- 
species Bose-Einstein condensate in a double well. Such 
a system is characterized by the following three factors. 
(1) The s-wave scattering — an and 022 for intra-species 
scattering and 0,12 for inter-species scattering; (2) Joseph- 
son tunnelling between the two wells; (3) population 
transfer between the two species driven by a resonant 
microwave field. We discuss the dynamics for three in- 
teresting regimes where we can obtain analytical results. 

(a) The strong coupling /tunnelling regime where the 
nonlinearity is negligible. We find the net effect of the 
coupling/tunnelling is either constructive or destructive; 

(b) the strong nonlinearity regime where Josephson tun- 
nelling is suppressed and the system behaves like a sim- 
ple two-level system; (c) the intermediate regime. For 
this case we only consider a specific example of vanish- 
ing intra-species scattering and weak inter-species scat- 
tering. We find collapses and revivals in the population 
dynamics. For the general case, we attack this problem 
numerically. We find the dynamics is rather complicated. 
The nonlinearity forces the system out of periodicity. We 
hope our work can be helpful to the study of quantum 
coherence of two-species BECs. 

We deeply appreciate Dr. F. Robicheaux at Auburn 
University and Dr. L. You at Georgia Institute of Tech- 
nology for enlightening discussions and helpful comments 
on the manuscript. We also thank J. Ludlow for proof 
reading our manuscript. This work is supported by 
grants with NSF. 
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